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ABSTRACT 

Context, cj Centauri is the most well studied Galactic Globular Cluster because of its numerous puzzling features: significant dis- 
persion in metallicity, multiple populations, triple main-sequence, horizontal branch morphology. He-rich population(s), and ex- 
tended star-formation history. Intensive spectroscopic follow-up observing campaigns targeting stars at different positions in the 
color-magnitude diagram promises to clarify some of these peculiarities. 

Aims. To be able to target cluster members reliably during spectroscopic surveys and both spatial and radial distributions in the cluster 
outskirts without including field stars, a high quality proper-motion catalog of a> Cen and membership probability determination are 
required. The only available wide field proper-motion catalog of ai Cen is derived from photographic plates, and only for stars brighter 
than 6 ~ 16. Using ESO archive data, we create a new, CCD-based, proper-motion catalog for this cluster, extending to 6 ~ 20. 
Methods. We used high precision astrometric software developed specifically for data acquired by WFI@2.2;;? telescope and pre- 
sented in the first paper of this series. We demonstrated previously that a 7 mas astrometric precision level can be achieved with this 
telescope and camera for well exposed stars in a single exposure, assuming an empirical PSF and a local transformation approach in 
measuring star displacements. 

Results. We achieved a good cluster-field separation with a temporal base-line of only four years. We corrected our photometry for 
sky-concentration effects. We provide calibrated photometry for UBVRcIc wide-band data plus narrow-band filter data centered on 
Ha for almost 360 000 stars. We confirm that the co Cen metal-poor and metal-rich components have the same proper motion, and 
demonstrate that the metal-intermediate component in addition exhibits the same mean motion as the other RGB stars. We provide 
membership probability determinations for published lo Cen variable star catalogs. 

Conclusions. Our catalog extends the proper-motion measurements to fainter than the cluster turn-off luminosity, and covers a wide 
area (~ 33' x 33') around the center of lo Cen. Our catalog is electronically available to the astronomical community. 

Key words. Globular clusters: general - Globular clusters: individual (to Cen) - Stars: populations II, H-R diagram - Catalogs - 
Astrometry 



1. Introduction 

The globular cluster co Centauri (w Cen) is the most luminous 
and massive cluster in the Galaxy. Observational evidence col- 
lected over the years has indicated that o) Cen is also the most 
puzzling stellar system in terms of stellar content, structure, and 
kinematics. Probably the most well studied of its peculiarities 
is one related to its stellar metallicity distribution (Norris & 
Besse ll [T975l [19771 Freem an & Rodgers [T975[ Bess ell & Norris 
[19761 Butle r et al. [19781 Norris & Da Costa [19951 Suntzeff & 
Kraft 119961 Norris et al. |1996l l. There is a significant disper- 
sion in the iron abundance distribution of co Cen, with a pri- 
mary peak about [Fe/H] ~ -1.7 — 1.8 and a long tail, extend- 
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ing to [Fe/H] ~ -0.6, which contains another 3-4 secondary 
peaks. It is possible to identify these metallicity peaks with dis- 
tinct stellar populations (Pancino et al. 120001 Rey et al. 120041 
Solli ma et a l. [20051 Villan ova et al. [20071 ). Ground-based (Lee 
et al. 119991 Pancino et al. 120001 ) and Hubble Space Telescope 
(HST) (Anderson p^T, Bedin et al. [20041 Ferraro et al. [20D4l ) 
photometry show clearly that o) Cen hosts different stellar pop- 
ulations. In particular, Pancino et al. ( I2000I I demonstrated that 
the CO Cen red giant branch (RGB) consists of at least four dis- 
tinct branches, spanning a wide range of metallicity. On the other 
hand, the oj Cen sub giant branch (SGB) has an intricate web of 
5 distinct sequences, indicating an extended range of metallicity 
and a ge (see Bedin et al. 120041 Hilker et a l. 120041 Sollima et al. 
[20051 Stanford et al.[2006l Villanov a et al . [2007l i . 

An derson (|1997| l. Bedin et al. (|2004| l. and Villanova et al. 
(120071 ). by studying fainter stars with deep and high resolution 
HST photometry, demonstrated that the main sequence (MS) is 
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Total WFI@2.2ra coverage around oiCen (1999 2003) 
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Fig. 1. Position footprint of the entire sample of WFI images 
around the lo Cen center (marked with a "+")• The first two rows 
show the covered areas sorted by month and year In the lower 
panels of the figure, the total coverage of all of the 279 images 
(on the left), and a zoom-in of the central part of the cluster 
(on the right) are shown. The numbers in parenthesis after the 
dates represent the total number of images for that observing 
run. North is up. East to the left. 



divided into 3 distinct sequences. The spectroscopic study of the 
MS stars of u Cen by Piotto et al. (2005) showed that the bluest 
of the main sequences is more metal rich than any of the redder 
sequences, which increased the ambiguity surrounding the clus- 
ter. An overabundance of He in the blue MS could reproduce the 
tj Cen MS photometric and spectroscopic properties (Bedin et 
al. 120041 NorrisEOQl Piotto et al. l2005T l. although the origin of 
the puzzling MS morphology is still far from being understood. 

A deeper insight into the enigmatic stellar populations of 
CO Cen should combine a deep, high resolution analysis of the 
inner and most crowded regions, with a wide field observations 
of the outskirts. While the first type of data has been provided ad- 
equately by HST, wide field coverage of o) Cen requires ground- 
based data which are more difficult to obtain. Acquisition of 
data for a wider field-of-view will inevitably result in higher 
contamination by Galactic foreground/background populations. 
The only reasonable and efficient way to decontaminate the 
to Cen outer stellar populations is by means of proper-motion 
analyses that help to isolate the Galactic contribution. The only 
available wide field co Cen proper-motion catalog (van Leeuwen 
et al. 120001 hereafter vLOO) is based on photographic observa- 
tions and only provides measurements for stars brighter than 
B ~ 16. In this paper, we attempt to provide the first CCD-based 
proper motion catalog of o) Cen, extending the cleaned stellar 
populations down to B ~ 20, i.e. 4 mag. deeper than vLOO. 

In the first paper of this series, Anderson et al. (120061) . here- 
after Paper I, demonstrated that WFI@2.2m observations, with 
a time base-line of only a few years, allow a successful separa- 
tion of cluster members from Galactic field stars in the two GCs 
closest to the Sun: NGC6121 and NGC6397. In this paper, we 



Table 1. Description of the data set used for the WFI@2.2m cat- 
alog. 



filter 


te.p 


seeing 


airmass 


January, 1999 


^842 


3x30s,lx300s; 


1".0-1".3 


~ 1.20 


658nm 


Ixl20s,5xl80s,lx900s; 


1".0-1".3 


1.15-1.18 




Mai-ch, 1999 






V843 


52x200s; 


0".8-l".2 


1.1-1.2 


^845 


51xl50s; 


0".7-l".5 


1.1-1.45 




April, 1999 








Ix5s,lxl0s,lxl5,lx30,5x60s; 


1".0-1".3 


1.3-1.6 


^845 


1 x5s, 1 X 10s, 1 x20s, 1 x45s,4x90s. 


0".74-l".7 


1.4-1.9 


658nm 


2x30,4xl20,5xl80s; 


0". 8-1". 15 


1.1-1.2 




(epoch I) 






g 


Ixl5s,lx30s,lx60s,5xl20s; 


0".75-l".3 


1.2-1.5 




3x5s,3xlOs,lxl5s,2x20s; 


0".7-l".0 


1.2-1.6 




Ix30s,4x45s,10x90s; 


0".76-l".36 


1.1-1.5 




July, 1999 






{J g^y 


2x 1800s; 


1".4-1".6 


~ 1.14 


-^842 


Ixl0s,lx30s,lx40s,lx300s; 


1".4-1".8 


1.14-1.25 


V843 


Ixl0s,lx20s,lxl50s,lx240s; 


1".3-1".6 


1.15-1.25 


-^845 


Ixl0s,2x20s,lx240s; 


1".13-1".6 


1.17-1.25 


658nm 


3x 120s, Ix 1200s; 


1".13-1".5 


1.18-1.23 




February, 2000 






^842 


2x30s,2x240s; 


1".45-1".7 


~ 1.13 


^843 


Ix30s,2x240s; 


1".1-1".2 


~ 1.13 


^845 


Ix30s,2x240s; 


~ 1".0 


~ 1.13 




March, 2000 






Um 


4x300s,lx2400s; 


1".1-1".2 


1.17-1.38 




5x30s; 


0".9-l".4 


~ 1.45 




June, 2002 






Uin 


4x30s,7x300s; 


0".8-2".0 


1.14-1.18 


Bgig 


Ix5s,3x8s,9x60s; 


0".8-l".5 


1.13-1.16 


^843 


3x5s,3x40s,3x60s; 


0".75-2".0 


~ 1.13 


^45 


3x20s,3x40s; 


0".7-l".4 


~ 1.13 




April, 2003 (epoch 11) 




^878 


7x40s;7 x 120s; 


0".7-0".9 


1.14-1.16 


Vs43 


7x40s;7 x 120s; 


0".8-l".0 


1.13 



apply the high precision astrometric and photometric techniques 
developed by Paper I to all available WFI@2.2;7z archive data of 
0) Cen. 

In Sect. 2, we describe the available WFI observations of 
(D Cen, and the data sets that we used to derive proper motions. 
In Sect. 3, we discuss our photometric data reduction technique, 
the sky-concentration effect minimization, and the photometric 
calibration. In Sect. 4, we describe in detail how we treated the 
differential chromatic refraction (DCR) effects between the two 
epochs. Membership probability is discussed in Sect. 5, while in 
Sect. 6 we outline possible applications of our catalog. Finally, 
in Sect. 7, we summarize our results and describe the electronic 
catalog. 
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Table 2. Characteristics of the used filters (from WFI manual) 
/lc=central wavelength, FWHM=Full Width at Half Maximum, 
/lp=transmission peak wavelength, Tp=transmission percentage 
at peak level. (*) LWP means Long Wave Pass: in this case the 
cutoff limit is determined by the CCD quantum efficiency. 



Wide-band filters 



Name 


A, 


FWHM 








[nm] 


[nm] 


[nm] 


[%] 


f/877 (U/50) 


340.4 


73.2 


350.3 


82.35 


f/841 (U/38) 


363.7 


38.3 


362.5 


51.6 


5878 (B/123) 


451.1 


133.5 


502.5 


88.5 


5842 (B99) 


456.3 


99.0 


475.0 


91.2 


V843 (V/89) 


539.6 


89.4 


523.0 


87.0 


RgM (Rc) 


651.7 


162.2 


668.5 


93.9 


/845 (Ic/lwp) 


783.8 


LWP* 


1001.0 


97.6 




Narrow-band filters 




Name 


A, 


FWHM 




Tp 




[nm] 


[nm] 


[nm] 


[%] 


658nm(H„) 


658.8 


10.3 


504.0 


90.7 



2. Observations 

We used a collection of 279 archive images acquired between 
January 20, 1999 and April 14, 2003 at the ESO/MPI2.2m tele- 
scope at La Silla (Chile) equipped with the wide-field imager 
camera (WFI). A detailed log of observations is reported in 
Table [1] This camera, which consists of an array mosaic of 
4x2 chips, 2141 x 4128 pixels each, has a total field of view 
of 34' X 33', and a pixel scale of 0.238"/pixel. More details of 
the instrumental setup were given in Paper I. Images were ob- 
tained using U, B, V, Rc, Ic wide-band and 658nm (He) narrow- 
band filters, whose characteristics are summarized in Table|2] 

For the derivation of proper motions, we used only B and 
V images acquired in April 1999 and April 2003 (see Sect. |4] 
for further details of this choice). The total field-of-view cov- 
ered by the entire sample is indicated in Fig. [T] where axis co- 
ordinates are (Aa cos 6, A6), expressed in units of arcmin from 
the (L> Cen center (North is up. East is to the left). Concentric 
circles have diameters, if not specified, of 10' (inner circle) 
and 30' (outer circle), and are centered on the cluster center: 
a = 201°.69065, 6 = -47°.47855 (Van de Ven et al. |2006l l. 

The first eight plots show the covered areas sorted by month 
and year, while in the bottom part of Fig. [T]the total coverage 
of all of the 279 images (on the left), and a zoom of the central 
part of the cluster (on the right) are shown. In the catalog pre- 
sented in this work, the proper motion measurements are avail- 
able only within the field-of-view in common between the two 
epochs used (see Fig.[T]i. 



3. Photometry, Astrometry, and Calibration 

3.1. Photometric reduction 

For the reduction of the WFI@2.2;7z photometric data, we 
used the software img2xym_WFI, a modified version of 
img2xym_WFC.09x10 (Anderson & King |2003T l. which was 
written originally for HST images, adapted successfully to 
ground-based data, and described in detail in Paper I. We closely 
followed the prescription given in Paper I for the data reduction 
of WFI images. This includes standard operations with the pixel 
data, such as de-biasing, flat-fielding, and correction for cosmic 
rays hits. 



At the basis of the star position and flux measurements, there 
is the fitting of the empirical Point Spread Function (PSF). In our 
approach, the PSF is represented by a look-up table on a very fine 
grid. It is well known that the shape of the PSF changes with po- 
sition in WFI@2.2ot chips. This variability can be modeled by 
an array of PSFs across the chip. The img2xym_WFI softwai-e 
works in a fully-automated way to find appropriate stars to rep- 
resent the PSF adequately. For practical purposes, the number 
of PSF stars per chip can vary between 1 and 15, depending on 
the richness of the star-field. An iterative process is designed to 
work from the brightest to the faintest stars and find their precise 
position and instrumental flux. A reasonably bright star can be 
measured with a precision of ~ 0.03 pixel (~ 7 mas) on a single 
exposure. 

Another problem of the WFI@2.2m imager is a large geo- 
metric distortion in the focal plane that effectively changes the 
pixel scale across the field-of-view. There are different ways to 
map this geometric distortion. We adopted a 9 x 17 element 
look-up table of corrections for each chip, derived from mul- 
tiple, optimally-dithered observations of the Galactic bulge in 
Baade's Window (Paper I). This look-up table provides the most 
accurate characterization of geometrical distortions available for 
the WFI@2.2m. At any given location on the detector, a bilinear 
interpolation between the four closest grid points on the look- 
up table provides the corrections for the target point. The de- 
rived look-up table may have a lower accuracy at the edges of 
a field, because of the way in which the self-calibration frames 
were dithered (see Paper I). An additional source of uncertainty 
is related to a possible instability distortions in the WFI@2.2m 
reported earlier. This prompted us to use the local transformation 
method to derive proper motions (see Sect.©. 

3.2. Si<y-concentration correction 

Once we obtained star positions and instrumental fluxes for all 
images, we had to minimize the so-called "sky-concentration" 
effect. The WFI@2.2m camera is affected significantly by this 
kind of light contamination (Manfroid & Selman 120011) . which 
is caused by spurious reflections of light at discontinuities in 
the optics and the subsequent redistribution of light in the focal 
plane. The insidiousness of the effect is due to the fact that this 
redistribution of light affects both the science and the flat-field 
exposures. 

Star fluxes are calculated by considering a local sky value, 
and therefore may be a negligible effect. However, since sky con- 
tamination also affects flat images, if it is not corrected properly 
during pre-reduction procedures, the quantum efficiency of the 
central pixel will be artificially lower with respect to that of the 
corner pixel. Consequently, the luminosity of a star measured 
in the middle of the mosaic camera will be underestimated by 
~ 0.1-0.2 magnitudes (in V band) with respect to the luminosity 
of the same star detected close to the mosaic edges. 

In Fig. 12^, we plot an instrumental color-magnitude diagram 
(CMD), which has been zoomed into the horizontal-branch (HB) 
region of co Cen, obtained by combining all V and Bms images 
of April 2003. We chose this particular data set to highlight the 
effect of sky concentration on un dithered images. In fact, this 
data set has only one pointing, as shown in Fig.[T] The positions 
of stars on the CCD mosaic is almost identical from one image to 



' Instrumental magnitudes are calculated to be -2.5 x log(S,DA',), 
where DN are the pixel's Digital Numbers above the local sky summed 
within a 10-pixel circular aperture. For a mean seeing of 0.8", saturation 
initiates at about 14.5. 
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Fig. 3. Final Amag correction grids for the different filters (B means the 8^42 filter only). Coordinates are in units of WFI pixels. Each 
WFI chip is highlighted by continuous black lines. Each element of the correction grids is colored according to the corresponding 
Amag correction applied. The grey scales are relative to the minimum/maximum correction for each filter. 



another, implying a small contribution to the r.m.s of the single 
star magnitude measurement due to sky-concentration effects. 
In Fig. 12^, we plotted 1252 stars with cry < 0.03 mag and cr^ < 
0.03 mag, where cry and erg are the standard errors of a single 
measurement (r.m.s). 

However, with only one pointing, sky concentration maxi- 
mizes its effect on the relative photometry of stars located at dif- 
ferent positions on the image. In Fig.|2^, we highlight the CMD 
of stars located at different positions on the CCD mosaic: with 
black dots (994 objects precisely), we show all stars between 
2' to 6' from a> Cen center, which is close to the mosaic cen- 
ter {x = 4150.69, y = 4049.97 on our master meta-chip). Red 
points are stars (258 objects) outside 12'. The displacement of 
the two HBs clearly shows that sky-concentration effects affect 
WFI photometry significantly if only one pointing is analyzed, 
and therefore needs to be corrected. 

In our case, the analyzed data sets for different filters come 
from several pointings (except the case of Z?878 images). In the 
process of matching all catalogs (for a given filter) to create a 
single master-frame, the true sky-concentration effect is reduced. 
For a given star, we considered the mean of the star fluxes orig- 
inating in different positions and for different pointings, so that 
the sky-concentration effect in the master frame was reduced. 



This process does, however, create systematic errors that affect 
the global photometry (see Selman l2001l l. 

In Fig. I2J3, we show the same zoomed HB region of co Cen 
derived, in this case, by matching all the available V and Z?842 
images obtained from the ESO archive. This data set contains 
several different pointings for both filters, so we were able to 
obtain photometry for the same stars located in some cases close 
to the mosaic center and in other pointings close to the mosaic 
edges. All plotted stars have again crBy < 0.03 mag. In this 
case, only 972 stars (with the same previous convention, 896 
black and 76 red) passed the selection criteria on the basis of 
photometric error. As explained before, matching catalogs for 
different pointing tends to minimize sky-concentration effects, 
but without an appropriate correction, rm.s. of measurements 
for the stars are enhanced. 



Andersen et al. (1 19951 1 studied the sky-concentration effect, 
typical of focal reducers, both by using simulations and analyz- 
ing data from the Danish telescope at La Silla. Their method 
for deriving the sky-concentration correction was based on the 
complex analysis of many star-field images taken at different ori- 
entations and positions during the night. Manfroid et al. (1200 11 1 
applied a similar method to derive the sky-concentration effect 
for WFI@2.2m data, while Selman ( 1200 11 1 developed a method 
to estimate the sky-concentration effect by the analysis of the ze- 
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V vs. {B,,,-V) [mstr,] 
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• 



ORIGINAL, ONE POINTING (2003) 

-0.4 -0.3 -0.2 



V VS. {Bg^r,-V) [instr.; 

0.3 -0.3 -0,1 

(b) • . 




ORIGINAL, ALL POINTINGS 

-0.3 -0.3 -0.1 




CORRECTED, ONE POINTING (3003) 

V vs. V) [instr.] 



CORRECTED, Al.l, POINTINGS 

V VS. {n„,^ V) [instr. 



Fig. 2. CMDs zoomed into the HB region of w Cen. All plot- 
ted stars have crsy < 0.03 mag. On the left, we used only 
images taken in April 2003 (i.e. only one pointing). On the 
right, we show photometry from all V and all ^842 images that 
have independent pointings. Upper diagrams are for the origi- 
nal photometric catalogs. The lower ones are derived from sky- 
concentration-corrected catalogs. Black dots are stars located 
close to the mosaic center, while red dots are stars close to the 
mosaic edges. 



ropoint variation in 3 dithered stellar frames, by evaluating this 
variations using a Chebyshev polynomial fit. They were able to 
reduce the internal error from 0.034 to 0.009 magnitudes in the 
V filter, and used the same polynomial fit to correct in addition 
the photometry in the other filters (see also Selman & Melnick 
12005b . Koch et al. (12004b provided an analogous prescription to 
correct for the sky-concentration effect by comparing photome- 
try derived with WFI and Sloan Digital Sky Survey (SDSS) data. 
Finally, Calamida et al. ( 12008b used some of the co Cen images 
that we present in this work to correct for the positional effects 
of the WFI camera by means of photometric comparisons with 
the local standard stars of ai Cen. 

The correction given by the ESO team, based on the V filter, 
consists of a 9* order bidimensional Chebyshev polynomial that 
should in principle be used also for the U and B filters. Selman 
( 12001b found that his solution for the V band was able to reduce 
the internal photometric error from 0.029 to 0.010 mag in B, and 
from 0.040 to 0.014 in U, while for the other filters the V cor- 
rection failed to reduce the internal photometric error Selman 
( 12001b argued that this is probably due to problems associated 
with atmospheric variations affecting data for different filters. 
Unfortunately, by using the same polynomial coefficients to cor- 
rect both V and B magnitudes, it is impossible to remove the 
color degeneracy due to the different response of the CCD to the 
sky concentration in the two different photometric bands. This 
degeneracy is of the order of ~ 0.04-0.05 mag in color in our V 
versus B - V CMD between inner stars (r ~ 4') and outer stars 
(r > 12'). 



Our adopted solution consists of a self-consistent autocali- 
bration of the sky-concentration map, and takes advantage of the 
high number of images analyzed, taken with different pointings. 
Below, we provide a description of the autocalibration proce- 
dure. We measured the raw magnitude mag, ^ of each /-star, in 
each y-image. We selected an image to be a reference frame (at 
the center of the dither pattern), and by using common stars with 
frame j we were able to compute the average magnitude shifts 
to bring each image-catalog onto the magnitude reference frame 
(A,). 

If there were no systematic errors, the same stars measured at 
different positions in two different frames, should have the same 
magnitude value, within the random measurement errors: 

mag,-o - Ao = . . . = mag,.^- - A^' = . . . = mag, ,,, - A„„ 

where m is the total number of images in that filter, used to per- 
form the autocalibration procedure. 

However, the same star closer to the center of the camera is 
systematically fainter than when it is measured closer to the cam- 
era edges. For each star, measured in several different frames, 
and at different positions on the camera, we can compute an av- 
erage of the values of the magnitudes in the reference system: 



/ 1 \ "' 



In the same way, we can compute a residual for the /-star in the 
;-image: 



(mag,-j - A;)-mag,.. 



All the residuals of stars close to the center will be systemat- 
ically positive, and those close to the edges systematically nega- 
tive. It appears clear that -at any given location on the camera- 
the average of the residuals from all the stars measured close to 
that location will provide a first spatial correction to our pho- 
tometry. It also appears clear that the determination of the sky- 
concentration photometric correction will be an iterative pro- 
cess. 

To guarantee convergence, we applied half of the recom- 
mended correction at the given location, to all our image cat- 
alogs. We then recomputed the Aj, and repeated the procedure 
until all the residual averages, at any given location, became 
smaller than 0.001 magnitudes. The null hypothesis of this pro- 
cedure is that the same star is imaged several times at different 
locations on the detector. To avoid systematic error, we also se- 
lect, for each filter, the same number of exposures per different 
pointing (as much as possible with the existing database). 

We use only those stars of the image catalogs with high S/N 
ratio (i.e. instrumental magnitude from -11 to -14) with a qual- 
ity PSF-fit smaller than 0.1 (as defined by Anderson et al. 12008b . 
and not too close to the cluster center (which does not necessar- 
ily coincide with the center of the camera) to avoid crowding, 
which compromises the photometric precision. The exact close- 
ness to the center of the cluster depends on the image exposure 
time; we excluded stars within a radius of 1000-2000 pixels from 
the cl) Cen center. 

At this point, we define the expression of "given location". 
After several tests, we found a spatial grid of 16 x 16 boxes to be 
the most suitable compromise between a large number of residu- 
als (dij) and a spatial resolution of correction sufficiently high to 
be useful. In Fig. [3] we show our final correction grids, respec- 
tively for U, ^842, V, Rc, Ic^ ™d Ha filters. For each filter, we 
present the final 16x16 element correction grid. Each element 
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Fig. 4. These figures show our calibration fits with the adopted color equations. For the B, V, Rc, Ic filters we used the on-line set of 
standard stars provided by Stetson, while for U and //„ bands we used as reference stars the Momany et al. (12003! ) catalog and the 
ACS/WFC catalog of Villanova et al. (i2007i i. respectively. See text for more details. 



is colored according to the grey scale values (black for the min- 
imum, white for the maximum). The grey scales vary linearly 
from the minimum to the maximum grid value for each filter It 
is clear that sky concentration affects different filters in different 
ways, and each filter must therefore be corrected independently. 
To evaluate the correction at any point of the camera, we com- 
pleted a bilinear interpolation of the closest 4 grid points. The 
adopted correction was less accurate close to the mosaic edges, 
where the peripheral grid-points have been stretched toward the 
boundaries. The available pointings for the other filters are also 
lower than for Bg42 and V, implying a less effective correction of 
the sky concentration. 

We emphasize that our correction is not completed by a star- 
to-star comparison. After the spatial correction, a single star can- 
not have a lower random error (higher precision) than mag,-. 
However, the systematic errors (accuracy) of single stars relies 
on the quality of individual grid-point solutions, which were al- 
ways calculated to be the average of residuals for several stars 
within each of the 16x16 cells. Even if our random errors for in- 



dividual stars are ~ 0.1 mag, with just 10 stars [in the worst case 
we still have at least 10 such stars] we can reduce our systematic 
errors to ~ 0.03 mag. The condition that each star has to be ob- 
served in each of the 16 x 16 cells is the ideal case. Deviations 
from this ideal case occur frequently, although overall, we are 
close to achieving the optimal solution. 

The total amplitude of our correction for the V filter is 0.13 
mag. For the same filter, Manfroid & Selman ( 12001b evaluated 
a 0.13 mag correction of similar spatial shape. Nevertheless, the 
two totally independent calibrations appear to be qualitatively 
the same. Based on the cell-to-cell scatter with the knowledge 
that sky-concentration is relatively flat, in V and / filters (for 
which we have more images) we estimate that the accuracy of 
our solution is as good as ~ 0.03 mag. Although our corrections 
do not use any color information, we note that the post-corrected 
CMD is in excellent agreement with (to within a few hundreds 
of a magnitude) the HB location (Fig.|2|. 

To verify qualitatively the high quality of our sky- 
concentration correction procedure, we show in Fig.|2]l the same 
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Fig. 5. {Left panels): plot of magnitude r.m.s (cr) as a function of magnitude. (Right panels): plot of the standard error of the mean 
{a I yjN - 1, where is the number of measurements) as a function of magnitude. 



CMD region (V vs. Z?842 - V), derived using all the available im- 
ages, after applying our correction. Stars located at different po- 
sitions on the meta-chip are not affected by the sky-concentration 
effect, and are located in the same CMD region. The total num- 
ber of plotted stars (1227, all with cr^ y < 0.03 mag) is compara- 
ble with that of Fig. |2^; this implies that we were able to remove 
the systematic contribution from our photometric rm.s values. 

Our solution works well for the available u Cen archive data 
sets (used to derive the corrections), but archive observations, 
in general, do not map every chip in a way that enables a sky- 
concentration correction that is universally applicable to be de- 
rived. We cannot guarantee that our solution can be applied to 
achieve the same positive results with other data sets. As proof 



of this issue, we applied our B^^o solution, to B^vg images with 
only one pointing. 

Images collected using the 8^42 and Bgyg filters are not 
so different, in term of central wavelength (see Table |2]): sky- 
concentration effects appear to be similar for almost identical 
filters (since they are related to atmospheric variations that af- 
fect the data for a range of different filters), so a photometric 
improvement is expected after correcting Bgvs images with our 
B842-derived solution. If our solution fails to correct the B%-j% 
photometry, this is probably due to the different pointings of 
the two B filters instead of the filters themselves. As shown in 
Fig. [2};, we found that a photometric improvement is present, 
with respect to Fig.|2^, but that the correction is not satisfactory. 
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Fig. 6. Photometric r.m.s for stars within 7' {top panel), from 7' 
to 14' {middle panel), and outside 14' {bottom panel). 



3.3. Instrumental UEVRdc-Ha photometric catalog 

We derived instrumental single-filter catalogs using all avail- 
able images, by matching each chip individually to minimize the 
zeropoint diff'erences between WFI chips. Included stars were 
measured in at least three distinct images. Photometric single- 
filter catalogs were then linked to the astrometric one. Linked 
star positions agreed with those in the proper-motion catalog 
within 1 pixel for BVRcIcHa filters, while for the U filter we 
had to adopt a larger matching radius (3.5 pixels). This is mainly 
due to the poorer distortion solution in the U band. 

Due to the aforementioned sky-concentration minimization 
problem with Bgvg images, our B photometry refers to the Bg42 
filter. 

As for the U photometry, we used only Usn images: in fact, 
Umi is a "medium" rather than a wide-band filter, of quite dif- 
ferent central wavelength and a low transmission efficiency with 
respect to t/gvv (Table |2]i. If not specified otherwise, we refer to 
Bg,42 and Uj^jj simply using B and U, respectively. 

3.4. Photometric calibration 

The photometric calibration of the WFI@2.2m data for BVRcIc 
bands was performed using a set of ~ 3000 online wide field 
photometric co Cen Secondary Standards stars (Stetson 120001 
12005b . The Secondary Standards star catalog covers an area of 
about 30' X 30' around the cluster center. We calibrated our U in- 
strumental photometry by cross-correlating our photometry with 
Momany et al. ( I2003I I U calibrated catalog (Stetson does not pro- 
vide U photometry for the o) Cen Secondary Standards). For Ha 
calibration, we used as reference-standard stars the 3x3 cen- 
tral ACS/WFC mosaic photometi-ic catalog in F658N band (GO 
9442), which was presented by Villanova et al. ( 2007 ). This HST 
catalog was obtained using img2xym_WFC.09x1 software; in- 
strumental magnitudes were transformed onto the ACS Vega- 



mag flight system following Bedin et al. (I2005I I. and by using 
the zero points of Sirianni et al. ( |2005t . 

For BVRcIc bands, we matched our instrumental magni- 
tudes and colors to the Stetson standard ones, and derived cal- 
ibration equations by means of an iterative least squares fitting 
of a straight line (see Fig.|4]i. For these filters, we found that only 
a first-order dependency of the color term affects our instrumen- 
tal magnitudes. The linearity of our calibration equations, which 
cover a wide range of colors (being derived from both HB and 
RGB stars) is evident from the plots of Fig.H] 

As in the calibration of our U instrumental photometry, 
the Momany et al. (120031) catalog was not corrected for sky- 
concentration effects. We found a magnitude dependence related 
to the star positions. We adopted a straight line fit to derive the 
calibration equation, because we were unable to consider the 
different color/magnitude dependencies individually. Therefore, 
our calibrated U magnitudes were not more reliable than the 0. 15 
magnitude level (maximal error). 

For the //„ filter, we again performed a straight line fit to 
derive the calibration equations, even though the data appear to 
suggest a second order color effect (see the corresponding panel 
of Fig.lHi. ACS/WFC data cover only the inner ~ lO'xlO' region 
of our catalog, and are therefore taken in extremely crowded 
conditions. This effect might strongly influence our photome- 
try mimicking the aforementioned second-order effect. In Fig.|5] 
we show in the left panels our photometric errors for each filter, 
as a function of the corresponding magnitude. The photometric 
errors (standard deviation) have been computed from multiple 
observations, all reduced to the common photometric reference 
frame in the chosen bandpass. In the right panels of Fig. |5] we 
plot the photometric standard error of the mean -to be defined as 
cr/ VA^ - 1, where is the total number of observations- versus 
the magnitude, for each filter 

To illustrate more clearly the dependence of our photomet- 
ric rm.s on crowding, we show (for the V filter only) in Fig|6] 
the photometric r.m.s o-y with respect to V for stars within 7' 
(top panel), from 7' to 14' (middle panel), and outside 14' (bot- 
tom panel). Stars located in the most crowded region of the field 
suffer higher uncertainty in their photometry. 

Our final catalog consists of about 360000 stars, in 
U BVRcIc wide-band and 658nm narrow -band filters, covering 
a wide area (~ 33' x 33') centered on o) Cen. We reach 3 magni- 
tudes in V band below the TO point with a photometric r.m.s of 
0.03 mag. 

3.5. Zeropoint residuals 

Even if our sky-concentration correction works well, residuals 
are still present, especially close to the corners of our final cat- 
alog. Due to the wide field area analyzed in this work, there is 
also the possibility of a contribution from differential redden- 
ing. Using ubvy Stromgren and V I photometry, Calamida et al. 
( 120051 ) developed an empirical method to estimate the differen- 
tial reddening of oj Cen. The authors found that the reddening 
can vary in the range 0.03 < E{B - V) < 0.15 from Stromgren 
filters, and 0.06 <E{B -V)<Q.l3 from V I filters, within their 
analyzed field-of-view of 14' x 14', which was centered on the 
center of co Cen. However, the results by Calamida et al. (120051 ) 
were questioned by Villanova et al. ( 120071 ). and the quantitative 
value of the differential reddening still needs to be confirmed. 

In the case of the B -V color, the maximum zeropoint resid- 
ual in our final catalog is less than 0.1 mag. To minimize any 
zeropoint variations, we used a method similar to that described 
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by Sarajedini et al. ( |2007| l. Briefly, we defined the fiducial ridge- 
line of the most metal-poor component of the o) Cen RGB and 
tabulate, at a grid of points across the field, how the observed 
stars in the vicinity of each grid point may lie systematically to 
the red or the blue of the fiducial sequence; this systematic color 
offset is indicative of the local differential reddening. 

Our online catalog magnitudes are not corrected for differ- 
ential reddening to enable the user to adopt their preferred cor- 
rection method in removing differential reddening and zeropoint 
residuals. 

4. Proper-motion measurements 

To complete the proper-motion analysis, we used only the B 
and V images taken in April 1999 (epoch I) and April 2003 
(epoch II). This choice was due to the fact that: (i) we have 
a fine-tuned geometric distortion correction map for V filter 
(Paper I), which has been proven to work well for the two B 
filters (Paper I); (ii) it offers the widest possible time base-line 
of ~ 4 yrs; and (iii) we have a relatively high number of images 
in both epochs, and with relatively deep exposures. 

We first photometrically selected probable cluster members 
in the V versus B - V color-magnitude diagram. These stars are 
located on the RGB (see the RGB selections in the top-panel of 
Fig. [13] within the magnitude interval 14.6 < V < 17.2. We 
used these stars only as a local reference frame to transform the 
coordinates from one image to the system of the other images at 
different epochs and therefore derive relative proper motions. By 
using predominantly cluster stars, we ensure that proper motions 
will be measured relative to the bulk motion of cluster stars. The 
expected intrinsic velocity dispersion of cj Cen stars for which 
we can measure reliable proper motions, is between 10 and 15 
km s"' (Merritt, Meylan & Mavor [T997| l. If we assume a distance 
of 5.5 kpc for oj Cen, as reported by Del Principe et al. (2006), 
and isotropic distribution of stars (good to first order), then these 
translate into an internal dispersion of 0.4-0.6 mas yr"'. 

Over the four-year epoch, the difference would result in a 
displacement of only 1.5-2.3 mas, which is a factor of 3 smaller 
than the random measurement errors (~ 7 mas). Conversely, the 
tangential velocity dispersion of field stars is a factor of ~ 10 
larger than the intrinsic velocity dispersion of the cluster. For 
field stars, proper motions are clearly not negligible with respect 
to measurement errors, and this has an adverse effect on the co- 
ordinate transformations. We removed iteratively stars from the 
preliminary photometric member list that had proper motions 
clearly inconsistent with cluster membership, even though their 
colors placed them close to the fiducial cluster sequence. 

To minimize the effects of geometric-distortion-solution 
residuals on proper motions, we used local transformations 
based on the closest 20 reference stars, typically extending over 
~ 30 arcsec. These were well-measured cluster stars of any mag- 
nitude selected to be on the same CCD chip, as long as their pre- 
liminary proper motion is consistent with cluster membership. 
No systematic errors larger than our random errors are visible 
close to the corners or edges of chips. 

To avoid possible filter-dependent systematic errors, we 
measured proper motions in the V and B bandpasses only, for 
which the geometrical distortion corrections were derived orig- 
inally (Paper I). Individual errors of proper motions for single 
stars were estimated as described in Sect. 7.3 of Paper I. For 
both epochs separately, we estimated the intra-epoch rm.s er- 
ror from all same-epoch plates transformed locally to the same 
reference frame. The proper-motion errors were computed to be 
the r.m.s of the proper motion, obtained by solving locally each 




Fig. 7. Proper-motion errors for stars within 7' {top panel), from 
7' to 14' (middle panel), and outside 14' (bottom panel). The 
proper-motion errors are expressed in the units of mas yr 



first-epoch frame into each second-epoch frame. These errors, 
however, were not entirely independent because the same frames 
were used more than once. Therefore, to obtain our most reli- 
able estimate of the proper-motion standard error, we added in 
quadrature the intra-epoch rm.s of each epoch. 

In Fig. |7] we show our proper-motion r.m.s, in units 
of masyr"', versus V magnitudes, calculated as cr^ = 



4 



cr^ cos (5 '^^'s■ '■°P P^iisl presents stars within 7' of the 
center of ca Cen, the middle panel is for stars between 7' and 14', 
while the lower panel shows the errors for stars outside 14'. The 
vertical dashed line indicates the saturation limit of the deepest 
exposures (V = 14.6), while the continuous line is at V = 16.5, 
the vLOO faintness limit. The precision of our proper-motion 
measurement is < 0.03 WFI pixels in 4 yrs down to V ~ 18 mag 
(i.e. cr < 1.9 mas yr"'). At fainter magnitudes, the errors gradu- 
ally increase, reaching ~ 5 masyr^'at V = 20. The stars brighter 
than y ~ 13 magnitude show a higher dispersion because of the 
image saturation even in the shortest exposures. Horizontal lines 
in Fig.|2]indicate the median proper-motion rm.s of unsaturated 
stars brighter than V = 16.5. We have 1.3 mas yr"', 1.1 masyr"', 
and 1.3 masyr 'for the top, middle, and bottom panel, respec- 
tively. The higher value for the inner stars is due to crowding 
while, for the outer stars, there is a combination of three factors: 
(1) our geometrical distortion solution is less accurate close to 
the WFI mosaic edges; (2) there are fewer cluster members, per 
unit area, usable as a reference for deriving proper motions; (3) 
we have a lower number of images that overlap with the external 
areas of the field-of-view. 

4.1. Cluster CMD decontamination 

To probe the effectiveness of our proper motions in separating 
cluster stars from the field stars, we show in Fig. [8] the vector- 
point diagrams (VPDs, top panels), and the CMDs in the V vs. 




Fig. 8. (Top-panels): proper-motion vector-point diagram. Zero point in VPD is the mean motion of cluster stars candidates. (Bottom 
panels): calibrated V, (B - V), color-magnitude diagram. (Left): the entire sample. (Center): stars in VPD with proper motion 
within 0.065 pixels (i.e. ~ 3.9 mas yr"') around the cluster mean. (Right): probable background/foreground field stars in the area of 
CO Cen studied in this paper. All plots show only stars with proper-motion cr smaller than 0.032 pixels (i.e. ~ 1.9 mas yr"') and V 
magnitude r.m.s. smaller than 0.02. 



B - V plane (bottom panels). In the left panels, we show the 
entire sample of stars; the middle panels display what we con- 
sidered to be probable cluster members; the right panels show 
predominantly the field stars. Plotted stars have a V r.m.s. lower 
than 0.03 mag. 

In the VPDs, we draw a circle around the cluster centroid of 
radius 3.9 mas yr '. Provisionally, we define as cluster members 
all points in the VPD within this circle. The chosen radius is the 
optimal compromise between missing cluster members with un- 
certain proper motions, and including field stars that have veloc- 



ities equal to the cluster mean proper motion. Even this approxi- 
mate separation between cluster and field stars demonstrates the 
power of proper motions derived in this study. A description of 
membership probabiUty is given in Sect.|5] 

4.2. Differential chromatic refraction (DCR) 

The DCR effect causes a shift in the photon positions in the 
CCD, which is proportional to their wavelength, and a function 
of the zenithal distance: blue photons will occupy a position that 
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Fig. 9. (Top panel): selected stars in the CMD of co Cen for DCR effect correction. (Middle and bottom panels): star displacements 
along X and Y axes, as a function of stars (B - V) color. The median shift of the nine samples are also showed, with errors. The 
continuous lines show the adopted fits used to quantify the DCR efl'ects. 



differs from that of red photons. The DCR effect is easier to de- 
tect and remove from CCDs, due to their linearity. 

Unfortunately, within each epoch, the available data sets are 
not optimized to perform the DCR correction directly (Monet et 
al. 119921 ). because the images have not been taken at indepen- 
dent zenithal distances. We can, however, check if possible dif- 
ferences in the DCR effect between the two epochs could gener- 
ate an apparent proper motion for blue stars relative to red stars. 

We selected four samples of stars located on the HB, and five 
on the RGB, as shown in Fig.|9](top panel), with different colors 
to estimate the DCR effect, in a magnitude interval of 1.8 mag 
in V (14.4<y<16.2), with proper motions < 3.8 masyr^'and 
r.m.s. < 1.9 masyr"'. We chose this magnitude range to: (i) 
avoid luminosity-dependent displacements (if any); (ii) include 



stars with a low r.m.s. in positions and fluxes (see Fig.|5]|7]i, ex- 
cluding saturated stars; and (iii) cover the widest possible color 
baseline according to the above points (i) and (ii). For each of 
the nine samples, we derived the median color and the median 
proper motion along ju„ cos 6 and ps, and their respective errors. 
Proper motions were expressed in the terms of a displacement 
over 4 years, that is in the units of WFI pixels along the X and Y 
axes of a detector (parallel to the RA and Dec directions). 

In Fig.|9] we show, in the top panel, the selected stars on the 
CMD used to examine the DCR effect; the linear fits adopted for 
the X and Y displacements are shown in the lower part of the 
figure. We found a negligible DCR effect along both X and Y 
axes. For this reason, we have not corrected our measurements 
for this effect. 
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Fig. 10. (Left): color-magnitude diagram, splitted into eight magnitude bins, for all the stars having a proper-motion r.m.s. increasing 
from 1.8 mas yr^'for the brightest bin to 5 mas yr"' for the faintest one, and photometric r.m.s. from 0.02 mag to 0.05 mag, generous 
enough to include main sequence stars down to V = 20. (Middle): proper motions for the stars in the corresponding magnitude 
intervals. A circle in each diagram shows the adopted membership criterion. (Right): color-magnitude diagram for assumed cluster 
members. 



In Fig.[TOl we show, on the left, our caUbrated V, B-V CMD, 
divided into eight magnitude bins. In each bin, we adopted dif- 
ferent selection criteria to identify cluster members, which were 
more stringent for stars with more reliable measurements from 
data of high signal-to-noise ratio, and less restrictive for star with 
less precise measurements. Plotted stars have a proper-motion 
r.m.s. of < 1.8, masyr"' for the brightest bin, to 5 masyr"' for 
the faintest one. The Photometric r.m.s. for both the bands range 
between 0.02 mag for the brightest bin to 0.05 for the faintest 



one, which is sufficient to include main-sequence stars down to 
V - 20. For each magnitude bin, we considered as cluster mem- 
bers those stars with a proper motion within the circle shown in 
the middle column of Fig. [10] 

On the right side of Fig. [TO] we show the color-magnitude 
diagram for stars assumed to be cluster members. The available 
archive images are again not sufficiently deep to derive reliable 
proper motions below the TO. The proper motions presented in 
this paper are not sufficiently accurate to study the internal mo- 
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tion of o) Cen. The main purpose of the proper motion presented 
in this work is to provide a reliable membership probability for 
spectroscopic follow-up projects, star counts, and the study of 
the radial distribution of the different branches (Bellini et al. in 
preparation). 

4.3. Astrometric calibration 

To translate the pixel coordinates into the equatorial coordinate 
system, we adopted the UCAC2 catalog (Zacharias et al. I2004I I 
as a reference frame. Due to the severe crowding in images of 
the inner parts of co Cen, this catalog was however inadequate 
for calibration purposes close to the center of ll) Cen (the central 
10' X 10' area corresponds almost entirely to a void in UCAC2). 
Another possible reference frame, especially for the cluster cen- 
ter, is the vLOO catalog. However, the precision of published co- 
ordinates is lower than ~ 20 mas and no analysis was provided 
by vLOO for the presence of potential systematic errors in the 
positions. Examination of vLOO proper motions by Platais et al. 
( |2003t indicated that a priori these systematic errors could not 
be discounted. These deficiencies in the vLOO positional cata- 
log were eliminated by re-reducing the original Cartesian coor- 
dinates (of formal precision equal to 2 mas), kindly provided by 
F. van Leeuwen. 

First, we selected only Class 0-1 stars from vLOO (i.e. their 
images were isolated or only slightly disturbed by an adja- 
cent image). Second, a trial equatorial solution was obtained 
for vLOO stars using the UCAC2 catalog as a reference frame. 
Third, the new set of vLOO coordinates was tested against the 
UCAC2 positions as a function of coordinates, magnitude, and 
color of stars. There are ~ 3000 stars in common between 
these two sets of coordinates. Assuming that the UCAC2 po- 
sitions are free of magnitude- color-related systematic errors, 
we found that the original vLOO Cartesian coordinates were bi- 
ased by up to 16 mas mag"' . There is also a detectable quadratic 
color-dependent bias along the declination. Both magnitude- and 
color-related biases were removed from the vLOO Cartesian co- 
ordinates before the final equatorial solution was obtained. 

The new reference catalog, covering a region of 1°5 x 1°5 
and magnitudes to V~16.5, contains 10291 stars and consists 
of approximately equal parts of the UCAC2 (trimmed down to 
stars with positional accuracies of higher quality than 75 mas) 
and the updated vLOO coordinates on the system of ICRS and 
epoch J2000.0. 

This new reference catalog was used to obtain the equato- 
rial coordinates of our co Cen stars. The WFl pixel coordinates 
of these stars were translated into global Cartesian system coor- 
dinates and corrected for geometric distortions. A simple low- 
term-dominated plate model was sufficient to calculate equa- 
torial coordinates. The standard error of this solution, employ- 
ing ~ 5500 reference stars, was 45-50 mas in each coordinate. 
These errors were higher than those listed in Yadav et al. ( 120081 
Paper 11), which is based on similar WF1@2.2ot data for the 
open cluster M67. We understand that image crowding remains 
a dominant source of increased scatter in our solution for o) Cen. 
Although we removed all stars with obviously poor astrometry, 
even a close but not overlapping image might slightly distort the 
position of a star, especially in photographic plates. The J2000 
positions of all stars for the epoch 2003.29 are given in Table 6. 

The proper motions in this work have not been translated 
into the absolute values, because there are too few background 
galaxies suitable for defining an absolute reference frame. 
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Fig. 11. (Top panels): vector-point diagrams of common stars in 
our catalog (left) and those of vLOO (right), with respect to the 
mean cluster motion. [Bottom panels): right ascension (left) and 
declination (right) proper-motion residuals as a function of the 
V magnitudes. Horizontal dashed and point-dashed lines show 
the 3cr-clipped median of the proper- mo tion dispersion for the 
selected cluster members, for our measurements and for vLOO, 
respectively. 



We compare our results with the proper-motion catalog by 
van Leeuwen ( |2000l l. First, we considered the common, unsat- 
urated stars in our catalog (V > 14.6) to the vLOO faint limit 
[V ~ 16.5). The selected samples contained ~ 3400 stars. Since 
vLOO proper motions are given in an absolute reference system, 
we subtracted from the individual vLOO proper motions the ab- 
solute mean motion of the cluster provided by the same authors 
[(//„ cos 6, //5)=(-3.97, -4.38) mas yr"' ]. 

In the top panels of Fig. (TT] we show on the left the vector- 
point diagram from our measurements, while on the right we 
show the vLOO values. In both diagrams, a concentration of stars 
at (0,0) masyr"' corresponds to the cluster members, while a 
diffuse distribution of stars around (-3, 5) masyr"', consists of 
the field objects in the foreground and possibly the background 
of o) Cen. 

The size of the proper-motion dispersions of the cluster 
members reflect both internal motions in the cluster and ran- 
dom errors. However, the internal motions are expected to be 
negligible. Assuming a distance of 5.5 kpc, and the Meylan et 
al. (|1995t measurement of the dispersion in the transverse ve- 
locity of ~10 kms ' (in the outskirts of the cluster that we are 
probing) the expected dispersion in the proper motions would be 
~ 0.4 mas yr"' . Our estimated errors for the selected sample was 
0.74 mas yr"' for cos 6, and 0.77 mas yr"' for jj.s. Therefore, 
the internal proper motions should not aff'ect more than 10-15% 
of the observed dispersions. 
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To estimate the observed proper-motion dispersion in the two 
samples (this paper and vLOO), we adopted the 68.27-th per- 
centile of the distribution (cr) about the median (estimated it- 
eratively with a 3cr-clipping), and for each coordinate indepen- 
dently. Due to the significance differences between cluster and 
field object motion, this procedure allowed us to isolate a sub 
sample of members. 

Our results were: 



this work : 



vLOO : 



crifXcr COS S) - 0.76 mas yr ' 
cr(yua) - 0.78 mas yr"' 

cr(iUa cos 6) = 0.94 mas yr"' 
crii^s) = 0.83 masyr"' 



For the selected sample, it is clear that our distribution is 
tighter, rounder, and in good agreement with our estimate of the 
errors. Even if our proper motions originate in images represent- 
ing half the total number of plates used by vLOO, we note that 
we study more than 3 mag fainter in V, and use a time baseline 
that equals only ~ 1/12 of that used by vLOO. 

The above performed test could be a bit unfair versus the 
vLOO catalog, because we used our non-saturated stars only, 
which are the faintest in the vLOO catalog. vLOO demonstrated 
emphatically that not all stars are suitable for astrometric mea- 
surements. We therefore performed a second test in which we 
chose the stars in vLOO with the most reliable measurements (be- 
longing to class and 1 only) that were brighter than V = 16, and 
had a good probability of being a cluster member (P^(vLOO) > 
75%); we compared the measurements for these stars with those 
in catalog. 

The results were as follows: 



this work : 



vLOO : 



crifia COS S) - 0.68 mas yr ' 
o-ijig) - 0.69 mas yr"' 

crifXaCosS) = 0.71 masyr"' 
cr(yua) - 0.62 masyr"' 



These results illustrate the slightly higher precision of the 
vLOO catalog, but the dispersions n the measurements are, in 
any cases, comparable. The dispersion obtained with our cat- 
alog, which now also includes saturated star measurements, is 
more reliable than found with the first test. The explanation of 
this apparent paradox is that, by selecting the best proper-motion 
measured stars in the vLOO catalog (class and 1, among the 
most isolated stars), we also selected the most isolated stars in 
our catalog, making the PSF-wings fitting more effective. Once 
again, it appears clear the huge potential that wide field CCD 
imager will have in astrometry in the future. 

5. Membership probability 

In the vector-point diagrams of Figs. [8] and [TO] two distinct 
groups of stars are clearly detectable: the bulk of stars belong 
to (L) Cen, with no mean motion (fia cos((5) = jUg = 0.0 mas yr"'), 
and there is a secondary broad group, which corresponds to field 
stars. 

Vasilevskis et al. (|1958t were the first to formulate the 
proper-motion membership probability. This method was later 
developed by many authors (Sanders 119711 Zhao & He 119901 
Tian et al. 119981 Balaguer-Nunez et al. 119981 and references 
therein) for several open and globular clusters. To derive our 
membership probability, we followed a method based on proper 
motions described by Balaguer-Nunez et al. ( |1998t . 



V 




12 



14 



16 



18 



20 



-0.5 




0.5 1 

(H-V) 



1 .5 



Fig. 12. Membership probability f ^ versus V magnitude for stars 
with cTyfi < 0.05 mag and proper-motion rm.s < 5 mas yr"'. 
Horizontal line show the = 90% level. The resulting CMD 
for stars with f ^ > 90% is shown on the bottom panel. 



First of all, we constructed frequency functions for both clus- 
ter and field stars, <J>'^. and respectively, derived from the clus- 
ter and field star distribution in the VPD. For the ;* star, these 
functions were as follows: 



exp 
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and 



exp 



(- 



1 

2(T~F) 



2;r(l-r2)l/2(^2. +^2)1/2(^2^+^2)1/2 



where (//i,,yUy,) are the /* star proper motions along the X and 
Y axes, (ex,, ey,) represent the respective displacement rm.s., 
(yU i^,/", / ) and (yU v^ ,/i, ,) are the central points of the field and clus- 
ter star proper motion, (ctv^ , cr,,^ ) and (ct v^ , cr,, ) are the field and 
cluster star proper-motion intrinsic dispersion, and y is the cor- 
relation coefficient, calculated to be 



7 = 



The spatial distribution was ignored due to the relatively 
small size of our field (~ 33' x 33') with respect to the oj Cen ra- 
dial extent of r, ^ 57' (Harris et al. 1 19961 1. For our calcula- 
tions, we considered only stars with an rm.s. in proper motion 
< 1.9 masyr 'to define <t)^ and O^. The center of the proper- 
motion distribution in the VPD for cluster stars was found to be 
at jCc = 0.00 masyr"', andy^ = 0.00 masyr"', with rm.s. values 
of (Tx, - 0.83 mas yr"', and cr,,^ - 0.83 mas yr"'. For field stars, 
we have: xj - -3.57 masyr"', yj - 5.12 masyr"', cTxf - 5.06 
masyr"', and cr,,^ = 2.86 masyr"', respectively. 

The distribution function for all the stars can be computed as 
follows: 

(D = (n, . (D^) + (nf ■ O}) = -h (Dj; 

where «c and n / are the normalized number of stars for cluster 
and field (n^ + itf - 1). Therefore, for the /* star the resulting 
membership probability is 



0.(0 

(D(0 ■ 



In Fig. [121 we show in the upper panel the distribu- 
tion versus V magnitude. To include also faint stars and reach 
V - 20, plotted stars have proper-motion r.m.s. < 5 masyr"' 
and cTyg < 0.05 mag. The horizontal line marks the 90% prob- 
ability level. The lower panel contains the CMD for stars with 

> 90%. 

When calculating formal membership probabilities for 
o) Cen, we have a paradoxical situation in which the main con- 
cern is to assign a reasonable probability to field stars. This is 
because in our sample the number of cluster stars is signifi- 
cantly higher than the small number of field stars. In addition, 
proper-motion errors have a strong dependence on magnitude 
(Figs. [71 [Toll, which is not accounted for in our membership prob- 
ability Pfj calculation. We therefore, also used the so-called local 
sample method (e. g.. Paper II) for membership probability cal- 
culation. In this method, for each target star a sub-sample of stars 
was selected to reflect closely the properties of a target. This as- 
sures a smooth transition in the calculated as a function of 
the chosen parameter. For w Cen, the obvious choice of parame- 
ter was the mean error cr^ of the proper motions. Given the wide 
range of cr^, we chose to consider only stars with cr^<7 mas yr"'. 
Less accurate proper motions are marginally useful for member- 
ship studies. We note that cr^ is calculated in the same way as in 
Sect. 4 and Fig. [H 

For each target star, we then selected a star a subsample of 
3000 stars almost with identical proper-motion errors to that of a 
target star The trial calculations indicated that we cannot model 
the distribution of field stars with a Gaussian because the number 



of potential field stars in the vector-point diagram is extremely 
low in the vicinity of the cluster-star centroid (Fig. [TTT i. A rea- 
sonable alternative to a Gaussian is a flat distribution of field 
stars. 

In essence, the membership probability P^, is driven by the 
distribution of cluster stars. If the modulus of the proper-motion 
vector of a star exceeds 2.5-3cr^, the corresponding is less 
than 1%. We provide these alternative estimates of P^ (called 
Pfi{2), to distinguish this from the first mentioned method Pf^il)) 
for 120259 stars. For the majority of stars, both values of P^ 
are similar. However, there are a number of cases in which the 
Pft estimates for the two methods differ radically. In the case 
of high proper-motion errors for cluster stars, the local sample 
method clearly provides a more realistic membership probabil- 
ity. A closer inspection of these cases indicates indirectly a po- 
tential problem in calculating the proper motion. If the error of 
proper motion along one axis is several times larger than the er- 
ror along the other axis, or if this error is too high for a particular 
magnitude, the chances are that our proper-motion value is cor- 
rupted and, hence, its formal membership probability is mean- 
ingless. 

Unless specified otherwise, we mean Pp(l) determination 
when referring to P^. 



6. Applications 

Our (D Cen proper-motion catalog can be used for different pur- 
poses. The first application was the selection of the most prob- 
able CO Cen members for spectroscopic follow-up studies. In 
Villanova et al. ( 120071 ). we used the proper-motion catalog of the 
present paper to pre-select sub-giant branch stars. This helped 
us to avoid the Galactic field stars close to the TO level. All 
resultant radial velocities were close to the oj Cen mean value, 
which confirmed their membership. On the other hand, the high 
photometric quality and the availability of several filters cover- 
ing a wide area around the cluster, imply that this catalog is an 
excellent photometric reference frame on which to register pho- 
tometry from different telescopes and cameras (Bellini et al. in 
preparation). 

Besides these obvious applications, our co Cen proper- 
motion catalog also provides an observational constraint of the 
origin of the composite stellar populations in oj Cen. Our cat- 
alog provides the necessary wide-field coverage and photomet- 
ric accuracy to investigate the radial distribution of the different 
o) Cen sub-population from the center of the cluster to ~ 22' 
(in the corners). We will report on this analysis in a subsequent 
paper (Bellini et al. in preparation). 

6.1. The proper motion of the RGB sub-populations 



Ferraro et al. ( 120021 1 cross-correlated the WFI photometric cata- 
log of Pancino et al. ( 120001 1 with the vLOO photographic proper- 
motion catalog. Their goal was to investigate the nature of the 
anomalous metal-rich RGB of o) Cen, the so-called RGB-a. 
In particular, they investigated the presence of proper-motion 
mean differences between the o) Cen bulk population (metal- 
poor RGB, so-called RGB-MP) and the minor, but yet impor- 
tant, metal-rich RGB-a population (Pancino et al. 120001 1. Their 
Fig. 2 showed significant variation in the relative proper mo- 
tion of RGB-a stars with respect to o) Cen RGB-MP motion. In 
particular, they found that the RGB-a stars had a mean proper 
motion of l^yUtotl - 0.8 masyr"' that is offset from that of the 
RGB-MP population. Therefore, they concluded that the RGB-a 
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subpopulation must have had an independent origin with respect 
to the RGB-MP one. 

Unsurprisingly, the Ferraro et al. (I2002I I study triggered new 
interest in the w Cen proper motion, and Platais et al. ( I2003I I 
presented a detailed reanalysis of the vLOO catalog. Platais et 
al. (I2003I I concluded that the reported proper-motion offset be- 
tween the 0) Cen sub-populations could be attributed to instru- 
mental bias. However, Hughes et al. (120041 ) commented that 
there was no residual color term in the omega Cen proper mo- 
tions, and that these authors misinterpreted the observed offsets. 
Specifically, Hughes et al. ( 120041 ) asserted that the summary ef- 
fect of color terms (before the corrections) amounted to no more 
than 1 masyr"' mag"' m B - V, while the offset in the proper 
motions for the anomalous omega Cen stars reached 2 mas yr"' 
and did not have the same direction as the color term. Regardless 
of the reason for the reported offset in the vLOO catalog, the pres- 
ence of this offset was not confirmed by the new HST-based pre- 
liminary proper motions (Anderson l2003l) . 

Two spectroscopic studies completed by Pancino et al. 
( I2007I I and Johnson et al. ( 120081 ) indicated that there was no evi- 
dence for any of the RGB stellar populations to have an offset in 
the mean radial velocity, or a different radial velocity dispersion. 
This result applied also to the RGB-a sub-population. 

Our astrometric catalog provided an independent data set 
with which we could test the Ferraro et al. (120021) conclu- 
sions. We repeated the same analysis performed in Ferraro et 
al. ( |2002| ). First, we divided the <jj Cen RGB population into 
three subpopulations (see top panel of Fig. [13]): RGB-MP (cyan), 
RGB-a (red), and RGB-Mint (RGB stars between RGB-MP and 
RGB-a, green). All plotted stars had high membership probabil- 
ity {Pf, > 90%), and photometric errors ranging from 0.02 mag 
for bright stars to 0.05 mag for faintest ones, in both filters. 

The RGB samples include all stars with V < 17.2. MP and 
Mint RGB stars merge with each other at fainter magnitudes. 
Sollima et al. ( I2005I I and Villanova et al. ( 120071 ) clearly showed 
that the SGB region of co Cen includes five subpopulations (see 
again Fig.fTSll. Three of these merge into each other and form the 
RGB-MP and RGB-Mint samples. We emphasize that our popu- 
lation sub-division extends to fainter magnitudes than in Ferraro 
et al. ( 120021 ). and therefore provides a higher RGB sampling and 
more robust statistics. Our samples contain 5182 RGB-MP stars, 
3127 RGB-Mint stars, and 313 RGB-a stars. 

A closer look at the three RGB sub-population proper mo- 
tions is presented in the bottom panels of Fig. [13] RGB-MP on 
the left, RGB-Mint in the middle, and RGB-a on the right. The 
red crosses report the median value (estimated iteratively) of the 
proper motion of the three sub-samples. For the RGB-MP sam- 
ple, we have: 

< fia cos 6 >- (0.00 ± 0.01) mas yr"' 

< i^s >= (-0.02 + 0.01) masyr-i, 

for RGB -Mint stars: 

< fia COS 6 >- (0.03 + 0.01) mas yr"' 

< Us >= (-0.01 + 0.01) masyr-i, 

and for the RGB-a sample: 

< fia COS 6 >= (0.05 + 0.03) mas yr"' 

< >= (-0.01 ± 0.04) masyr '. 

We found no evidence for the presence of differences among 
the relative proper motions of the three RGB sub-populations at 
the level of 0.05 mas yr"' in yU„ cos 6 and of 0.04 mas yr"' in fig 
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Fig. 13. (Top panel): V vs. B-VCMD of proper-motion-selected 
member stars on which we defined three RGB subsamples: 
RGB-MP (cyan), RGB-Mint (green), and RGB-a (red). See text 
for the definition of the subpopulations. {Bottom panels): Proper 
motions of RGB-MP (left), RGB-Mint (middle), and RGB-a 
(right) subpopulations. 



(i.e. relative tangential velocities of < 1.3, and < 1.1 km s"', 
assuming a distance of 5.5 kpc for co Cen). All three RGB sub- 
samples exhibit the same mean proper motion, within the errors. 

We therefore, agree with the results of both Platais et al. 
(|2003T l and Johnson et al. ( |2008] l for RGB-a, and show that 
the RGB-Mint proper motion is also consistent with the other 
o) Cen sub-populations. 

A final word on this issue requires an internal stellar proper- 
motion analysis, but suitable catalogs are not yet available. 

6.2. Membership probability of published d Cen variable 
stars 

The study of variable stars in oj Cen can certainly benefit from 
our proper-motion catalog and cluster membership derivation. 

Using the 1 .Qm Swope Telescope, Kaluzny et al. ( 120041 ) mea- 
sured light curves of ~ 400 variable stars in oj Cen, 1 17 of which 
were new identifications. We cross-checked our proper-motion 
catalog with their own, and found a total of 338 variable stars in 
common, which had both V and B measurements in our catalog. 
In particular, there were 90 variable stars for which Kaluzny et 
al. ( 120041 ) did not provide B and/or V magnitudes. Our proper- 
motion catalog was also useful to locate these stars in the color- 
magnitude diagram. The position of the Kaluzny et al. ( 120041 ) 
variables in our CMD, as well as their membership probability, 
is shown in the left panel of Fig. [14] Depending on our measured 
membership probability, we divide the Kaluzny et al. (2004]) 
sample into three categories: (1) < 30% (red squares); (2) 
30% < P^, < 80% (green triangles); and (3) > 80% (blue 
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Fig. 14. {Left): to Cen CMD with the cross-checked variable stars from the Kaluzny et al. (I2004I I catalog; red squares mark stars 
with a membership probability < 30; green triangles are stars with 30 < < 80, while blue circles are those with a probability 
to be cluster members > 80. The membership probability versus the V magnitude is also shown. (Right): comparison with 
Weldrake et al. ( I2007I I catalog. 



circles). Of the 1 17 new variable identifications, 1 12 are in com- 
mon with our catalog. Of these 112 stars, 15 have < 30%, 
and therefore these are probably not cluster members. On the 
other hand, 19 of these new variables have 30% < P/j < 80% 
and their membership remains uncertain. The remaining 78 stars 
(Pf^ > 80%) are almost certainly oj Cen members. In Table[3] we 
report the membership probability values for all Kaluzny et al. 
(12004b variable stars identified in our catalog. New Kaluzny et 
al. ( 120041 ) identifications have their IDs starting with "N". 

The Kaluzny et al. ( 120041 ) variable star catalog was cross- 
checked with the X-ray sources detected by the XMM-Newton 
analysis presented by Gendre et al. ( 120031 ). For the 9 stars in 
common (see also Table 2 of Kaluzny et al. 120041 ). we can now 
provide (see Table HI a membership probability based on our 
proper-motion analysis. A quick glance at Table 4 allows us to 
infer that only NV383 and NV369 are very likely cluster mem- 
bers, whereas the remaining 7 stars are most probably field pop- 
ulation. 

Using the 1.0m Telescope of the Australian National 
University, Weldrake et al. (12007b detected a total of 187 vari- 
able stars covering a wide area around a> Cen. These stars were 
matched with the Kaluzny et al. ( 12004b catalog, and 81 vari- 
able stars were found to be new discoveries. A cross-check of 
our catalog with that of Weldrake et al. ( 12007b provided 102 



variable stars in common. Their location in the V versus V - I 
CMD is shown in the right panel of Fig. [14] As completed for 
the Kaluzny et al. (12004b sample, different symbols mark dif- 
ferent membership probability ranges, and Table |5] reports our 
membership probability for the Weldrake et al. ( 120071 ) variable 
catalog. Of the 81 Weldrake et al. ( 12007b new variable stars, 16 
have counterparts in our proper-motion catalog. Four of these 
new variables are clearly not cluster members. These field ob- 
jects include a detached eclipsing binary (V59), and three long- 
period variables (V31, V125, and V126). 



Table 4. Membership probability of XMM-Newton X-Ray 
counterparts candidates. The IDs refer to Kaluzny et al. work. 



ID,, 


P, 




P, 


NV367 





NV376 





NV375 





NV377 





V167 


23 


NV369 


94 


V223 





NV383 


99 


NV378 


31 
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Table 3. Membership probability for the Kaluzny et al. ( |2004| l variable star catalog. ID/c are the Kaluzny identification labels, while 
ID,„ refer to this work. 



IDk 




ID,,, 


IDk 




ID,,,, 






ID,,,. 


IDk 


P, 


ID,,,. 


IDk 


P, 


ID,„ 


IDk 


P, 


ID,,, 


VI 


99 


290466 


V66 


99 


311494 


V132 


96 


184072 


V214 


97 


42989 


NV296 


99 


356415 


NV354 


97 


263543 


V2 


74 


271987 


V67 


99 


352882 


V135 


97 


181199 


V216 


99 


155556 


NV297 


32 


250290 


NV355 


93 


196052 


V3 


99 


248511 


V68 


100 


345639 


V136 


99 


213549 


V217 


98 


218342 


NV298 


100 


87382 


NV356 


49 


206069 


V4 


97 


279540 


V70 


98 


93459 


V137 


99 


225149 


V218 


96 


358729 


NV299 


98 


310461 


NV357 


92 


158687 


V5 


100 


298983 


V71 


99 


209884 


V139 


90 


215311 


V219 


88 


347441 


NV300 


98 


284928 


NV358 


99 


129801 


V8 


99 


200194 


V74 


100 


360819 


V141 


100 


177819 


V220 


100 


321724 


NV301 


96 


170554 


NV359 


46 


222238 


V9 


100 


238167 


V75 


99 


352935 


V142 


65 


193092 


V221 


96 


320315 


NV302 


19 


252220 


NV360 


100 


118132 


VIO 


98 


273997 


V76 


100 


340408 


V143 


90 


217472 


V223 





355759 


NV303 


98 


145690 


NV361 


64 


234998 


Vll 


99 


302079 


V77 


99 


316247 


V144 


95 


201059 


V224 


99 


349215 


NV304 


98 


341374 


NV363 


93 


321474 


V12 


99 


283465 


V78 





243506 


V145 


91 


143665 


V225 


96 


274338 


NV305 





343091 


NV364 


93 


243121 


V13 


99 


259326 


V81 


99 


270509 


V146 


72 


177889 


V226 


99 


167657 


NV306 


96 


198996 


NV365 





71679 


V14 


97 


26311 


V82 


99 


236315 


V147 


97 


143269 


V227 


98 


287148 


NV307 


82 


211960 


NV366 


97 


132355 


V15 


100 


273754 


V83 


99 


323500 


V148 


95 


209005 


V228 


96 


92375 


NV308 


98 


327308 


NV367 





350053 


V16 


100 


40312 


V85 


98 


292312 


V150 


100 


58317 


V229 


97 


285982 


NV309 


85 


193994 


NV368 





349013 


V17 


99 


261135 


V86 


99 


242948 


V152 


82 


177559 


V231 


99 


59688 


NV310 


15 


127228 


NV369 


94 


124330 


V18 


98 


267824 


V87 


97 


255154 


V153 


100 


238719 


V233 


99 


336053 


NV311 


91 


170132 


NV370 





358442 


V19 


100 


205260 


V88 


94 


261450 


V154 


96 


155744 


V234 





157490 


NV312 


95 


183516 


NV372 


93 


290428 


V20 


99 


204861 


V90 


89 


239287 


V155 


98 


272492 


V235 


87 


309505 


NV313 


98 


134090 


NV373 


59 


300775 


V21 


99 


246994 


V91 


92 


241579 


V156 


99 


128871 


V236 


85 


282601 


NV314 


85 


163946 


NV374 


99 


312254 


V22 


99 


86072 


V92 


99 


327526 


V157 


98 


221081 


V237 


99 


286453 


NV315 


60 


267539 


NV375 





349760 


V23 


100 


273446 


V94 


98 


305927 


V158 


31 


153087 


V238 


99 


179489 


NV316 


99 


142842 


NV376 





189233 


V24 


99 


84169 


V95 


94 


188847 


V161 


76 


149737 


V239 


99 


179912 


NV317 


100 


361100 


NV377 





324933 


V25 


99 


199162 


V96 


63 


225607 


V162 


97 


108440 


V240 


99 


239985 


NV318 


98 


246801 


NV378 


31 


341736 


V26 


98 


226616 


V97 


99 


256623 


V163 


100 


337663 


V241 


1 


289203 


NV319 


94 


145544 


NV379 


99 


328605 


V27 


96 


201410 


V98 


99 


227653 


V164 


99 


334394 


V242 


98 


171250 


NV320 


99 


80363 


NV380 


96 


148727 


V29 


62 


191401 


V99 


99 


210796 


V165 


85 


227981 


V249 


99 


236550 


NV321 


87 


209026 


NV381 





311380 


V30 


99 


168094 


VI 00 


77 


215927 


V166 


72 


241584 


V250 


92 


62410 


NV322 


93 


214126 


NV382 


96 


343538 


V32 


100 


322448 


VlOl 


99 


169863 


V167 


23 


87332 


V251 


97 


231812 


NV323 


65 


265232 


NV383 


99 


336713 


V33 


98 


184952 


V102 


100 


162564 


V168 
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7. Summary 

We have applied the photometric and astrometric technique de- 
veloped by Anderson et al. (2006, Paper 1) to the most puzzling 
globular cluster of the Milky Way: to Centauri. 

Based on CCD observations taken with only four years 
of temporal base-line, our measurements provide accurate 
proper motions to B ~ 20, four magnitudes deeper than the 



photographic catalog of vLOO. We have minimized the sky- 
concentration effects in our photometry. We provide a mem- 
bership probability for all stars. Our catalog contains almost 
360 000 stars with measured proper motion, and covers a wide 
area (~ 33 x 33 arcmin-) around the cluster center. In Fig. [15] 
we show a summary of our photometric catalog; we plot several 
ai Cen CMDs, derived with all the available filters and different 
color-baselines (top and middle rows). Plotted stars have a mem- 
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Table 5. Membership probability for the Weldrake et al. ( 120071 ) variable star catalog. IDw are Weldrake et al. identification labels, 
ID,„, come from this work. The symbol (*) marks new Weldrake et al. identifications. 







ID,,,, 
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1D„ 




ID,,,, 
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96 
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352935 


VlOl 


99 


36832 
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35 


151268 
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51 


218172 
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93 


287079 
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36 
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bership probability of P^, > 90%. Photometric errors range from 
0.02 mag for brighter stars to 0.05 mag for the fainter ones. In the 
bottom panels of Fig. [15] we show the SGB region of ca Cen in 
the B versus B-R CMD on the left, and two color-color diagrams 
with different colors on the right. 

The high precision of our astrometry and multi-band pho- 
tometry once again emphasizes the importance of accurate rep- 
resentation of the PSF across the entire field-of view, especially 
for wide-field imagers, exemplified by the concept of empirical 
PSF (Paper I). 

The primary aim of this work is, of course, to provide wide- 
field membership probability measurements for spectroscopic 
follow-up studies, down to the turn-off region of the cluster 
However, the high quality of our photometric and astrometric 
measurements also provide a crucial observational constraint of 
the multiple ai Cen sub-populations. Due to our proper-motion- 
selected RGB sub-populations, we can confirm that the metal- 
poor, metal-intermediate, and metal-rich components have the 
same proper motions which is that of co Cen, within our mea- 
surements uncertainties. 

We finally provide membership probability determinations 
for the Kaluzny et al. ( |2004] l and Weldrake et al. ( |2007] l 
o) Cen variable star catalogs. 

7.1. Electronic catalog 

The catalog is available at the SIMBAD on-line databas^. 
Description of the catalog: column (1) contains the ID; columns 
(2) and (3) give the J2003.29 equatorial coordinates in decimal 
degrees; columns (4) and (5) provide the pixel coordinates x and 
y of the distortion-corrected reference meta-chip. Columns (6) 
and (7) contain the proper-motion values along /i„ cos S in units 
of masyr ', with the corresponding rm.s.; columns (8) and (9) 
provide the proper-motion values along /y^, with the correspond- 
ing r.m.s., in the same units. Columns from (10) to (21) give the 
photometric data, i.e. U, B,V,Rc, Ic, Ha magnitudes and their 
errors. The last two columns give the proper- motion member- 
ship probability P/jil) for all the stars (22), and for a sub-sample 
of 120,259 stars with the alternative membership determination 
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Pfi{2) (23). In this case, if the second membership determination 
is not provided, we flag it with the value -1. 

If photometry in a specific band is not available, a flag equal 
to 99.999 is set for magnitude and error. 

Table 6 contains the first lines of the on-line catalog. 
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20.116 


0.017 


19.692 


0.068 


19.347 


0.069 


99.999 


99.999 


62 




9 


201.714951 


-47.750153 


3904.018 


-91.634 


5.27 


3.52 


2.33 


6.89 


99.999 


99.99 


9 20.354 


0.050 


19.696 


0.028 


19.261 


0.051 


18.909 


0.032 


99.999 


99.999 


70 




10 
[...] 


201.629034 
[...] 


-47.750053 
[...] 


4779.417 
[...] 


-91.229 
[...] 


4.24 
[...] 


3.80 
[...] 


7.46 
[...] 


1.69 
[...] 


99.999 
[...] 


99.99 
[...] 


9 19.275 
[...] 


0.024 
[...] 


18.694 
[...] 


0.024 
[...] 


18.310 
[...] 


0.030 
[...] 


17.953 
[...] 


0.057 
[...] 


99.999 
[...] 


99.999 
[...] 



[...] 


55 
[...] 



Table 6. First lines of the electronically available catalog. 



